arXiv:1502.01964vl [cs.NI] 6 Feb 2015 


Maximum Likelihood based Multihop Localization 
in Wireless Sensor Networks 


CamLy Nguyen\ Orestis Georgiou^, and Yusuke Doi^ 

^Network System Laboratory, Corporate Research & Development Center, Toshiba Corporation, 
1 Komukai-Toshiba-cho, Saiwai-ku, Kawasaki 212-8582, Japan 
^Toshiba Telecommunications Research Laboratory, 32 Queens Square, Bristol, BSl 4ND, UK 


Abstract —For data sets retrieved from wireless sensors to be 
insightful, it is often of paramount importance that the data 
be accurate and also location stamped. This paper describes a 
maximum-likelihood based multihop localization algorithm called 
kHopLoc for use in wireless sensor networks that is strong in both 
isotropic and anisotropic network deployment regions. During 
an initial training phase, a Monte Carlo simulation is utilized 
to produce multihop connection density functions. Then, sensor 
node locations are estimated by maximizing local likelihood 
functions of the hop counts to anchor nodes. Compared to other 
multihop localization algorithms, the proposed kHopLoc algo¬ 
rithm achieves higher accuracy in varying network configurations 
and connection link-models. 

Index Terms —Localization, range-free, wireless sensor net¬ 
works, mesh networks, multihop, connectivity. 

1. Introduction 

Wireless sensor networks (WSN) are composed of a set 
of spatially distributed wireless nodes, with sensing and 
transceiving capabilities, tasked with monitoring physical or 
environmental conditions, such as temperature, sound, pres¬ 
sure, radiation, etc. Data collected is then wirelessly passed 
through the network to a main gateway for storing and 
processing. WSNs are essential for applications such as en¬ 
vironmental monitoring, target tracking, disaster relief and 
rescue operations 0,0. Nowadays, they are also becoming 
an indispensable part of smart technologies with applications 
in smart cities and smart buildings 0- 

In view of the internet of things future vision, almost every 
device will soon have transceiving capabilities, be packed with 
sensors, connected to a network and producing huge data sets. 
As is often the case however, location information is vital for 
the insightful processing of this data. GPS modules may be 
embedded to each sensor node enabling it to autonomously 
discover its location both accurately and on demand, however 
this does not come without a cost to the manufacturer (and 
hence the user) and the node’s power source - about 30mA 
at 3.3V. Moreover, in some extreme instances such as sand 
storms and blizzards, or simply when operating indoors, 
satellite signals cannot reach the sensor nodes. 

To alleviate such problems, cooperative schemes have been 
developed to estimate the locations of sensor nodes with the 
assistance of nodes which have perfect location information 
0, ©• The nodes whose locations are known and the nodes 
whose locations are unknown are usually called anchor nodes 


and target nodes respectively and the localization techniques 
can be broadly classified into two schemes: range-based 
and range-free schemes. Range-based schemes assume 
that the distance or angle between anchor nodes and target 
node can be measured based on signal measurements such 
as received signal strength indication (RSSI), time of arrival 
(TOA), or angle of arrival (AOA). In large-scale WSNs where 
signal range is limited however, range based schemes typically 
require a lot of anchor nodes to produce accurate results. On 
the other hand, range-free schemes 0 estimate inter-node 
distances based on hop count information, thus all target nodes 
can be localized with fewer anchor nodes. 

Conventional range-free approaches 0 -p0| usually con¬ 
sider isotropic WSNs where sensor nodes are uniformly dis¬ 
tributed in a regular region (e.g. a square domain), thus the 
distance between two nodes is assumed to be proportional 
to their hop count. The celebrated DV-hop algorithm esti¬ 
mates the average one-hop distance, and then multiplies this by 
the hop count to at least three anchor nodes before trilaterating. 
Although improved DV-hop algorithms have been suggested 
in the literature Q, pQ| , performance gains have been limited, 
even more so when used in anisotropic WSNs where factors 
such as irregular radio propagation, obstacles, nonuniform 
node distributions degrade the hop-distance proportionality as¬ 
sumption. DV-hop-like variants which are anisotropic network 
compatible have also been proposed (TT], fT^ , the main idea 
usually being to reduce the estimation error by reducing the 
effect of unreliable anchors. Thus, these algorithms require an 
increased number of anchor nodes, possibly as many as range- 
based schemes. Other approaches (T3), |T4) attempt to use 
approximate shortest paths to reduce the effect of anisotropic 
networks, yet jfMj under performs in irregular-shaped regions 
as shown in |13| which comes with large communication and 
computational overhead. 

In this paper, we propose a maximum likelihood based mul¬ 
tihop localization algorithm called kHopLoc which achieves 
good performance in both isotropic and anisotropic WSNs. 
The algorithm first runs a training phase during which a Monte 
Carlo simulation is utilized to produce accurate multihop 
connection probability density functions (described later). In 
its second phase, the algorithm constructs likelihood functions 
for each target node based on their hop counts to all reachable 
anchor nodes which it then maximizes to produce localization 


information. Unlike most DV-hop algorithms which use only 
first order statistics, the proposed kHopLoc algorithm gen¬ 
erates and uses the full multihop density distributions (even 
for anisotropic networks) thus constructing accurate likelihood 
functions and in turn localization results. In addition, our 
algorithm’s communication cost is about half of most DV- 
hop-like algorithms and computational cost is much smaller. 

The rest of this paper is organized as follows. In Section 
l^we describe the network and system model. In Sec. Ill 


we 


present the details of our proposed kHopLoc algorithm. In 
Sec. IIV] we evaluate the performance of kHopLoc through 
numerical simulations. Finally, Sec. [V| concludes the paper 
with a summary and some discussion on future work plans. 
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Fig. 1. Plots of the pair connectedness function H using (a) the Rayleigh 
fading communication model for parameters /3 = 1 and 77 = 2, 3, 4, 
and (b) the Quasi Unit Disk Graph communication model for parameters 
dmax = 1 and DOI = 2 , 3, 4. 


function can be expressed as 


II. Network Definitions and System Model 


H = exp{-^d^), ( 2 ) 


Consider a WSN of N total sensor nodes, M of which 
are anchor nodes (i.e. have perfect location information). The 
remaining N — M sensors are target nodes (i.e. locations 
are unknown). All nodes are randomly distributed in some 
subset of and are equipped with isotropic antennas. Due to 
phenomena such as fading and multi path, the communication 
model adopted may not be well described by the simplistic 
disk model (where two nodes connect if they are within a finite 
range of each other). Instead, we adopt a random connection 
model where nodes connect with a distance dependent proba¬ 
bility while also accommodating for environmental parameters 
such as path loss exponent GD- Namely, we consider two such 
communication models in order to show that the kHopLoc 
algorithm works in different communication models: Rayleigh 
fading communication model and Quasi Unit Disk Graph 
(QUDG) communication model GD- The latter model also 
allows the direct comparison of kHopLoc with that given 
in GD- Moreover, we consider isotropic networks in which 
the sensor nodes are distributed in a square, and anisotropic 
networks where sensor nodes are deployed in irregular shaped 
regions. 

A. Rayleigh fading communication model 

The pair connectedness function H defined as the proba¬ 
bility that two nodes are directly connected. One way of for¬ 
mulating H is thus the complement of the information outage 
probability with respect to a mutual information threshold d. 
For a narrow band transmission subject to small-scale fading 

ff = P(log2(l + SNR-X)>^)=Fx(^^^), (1) 

where X is a random variable defining the gain of the wireless 
channel between the two nodes, Fx is its complementary 
cumulative distribution function (CCDF), and SNR is the long¬ 
term average received signal-to-noise ratio. Since the nodes are 
equipped with isotropic antennas, the SNR oc d~^, where d 
is the distance between the two nodes, and p is the path loss 
exponent. Adopting a Rayleigh fading model for the small- 
scale fading gain X ^ exp(l) implies that the CCDF of X 
is written as Fx{x) = exp(—x), and the pair connectedness 


where /3 is a constant depending on the node transmission 
power, wavelength, information threshold etc. and can be 
understood as an effective communication range vq = 

Notice that in the theoretical limit of 77 ^ 00 , we have that 
ro ^ 1 and H converges to the unit disk model. Without loss 
of generality, we will henceforth set /3 = 1. 


B. Quasi Unit Disk Graph communication model 

In the QUDC communication model (T^, the probability 
H that two nodes a distance d apart are directly connected is 


/ 

1 

DOI{dmax-d) 

dmaxiDOI-1) 

0 


if d< dmax/DOI 

if d G [dmax / DO I ^ dmax] 

if d> dmax/DOI 


(3) 


where dmax > 0 and DOI > 1 are the maximum successful 
transmission distance and the degree of radio irregularity 
respectively. 

During computer simulations, the network graph edges are 
formed if a random number ( G [ 0 , 1 ] is less than the 
calculated H of the respective communication model. Two 
nodes sharing a successful link are called one-hop neighbours. 
More generally, two nodes with the minimum hop count of 
k > 1 (measured along the shortest path) are called k-hop 
neighbours. The two H functions are plotted in Fig. 


HI. kHopLoc algorithm description 

Our proposed kHopLoc algorithm is composed of three 
simple steps. In the first step, each target node counts the 
minimum number of hops to each anchor node. In the second 
step, the conditional probability density function is generated 
for the probability density function p{d\k) of the Euclidean 
inter-node distance d given a hop count from target to anchor 
node equal to /c = 1, 2,.... This step is preferably performed 
centrally using a Monte Carlo simulation after which fitting 
parameters are fiooded through the network. Finally, in the 
third step, each target node calculates its own position by max¬ 
imising the joint conditional density function to all reachable 
anchors. We now describe each of these steps in more detail. 
















A. Step 1: Minimum hop-count to anchor nodes 

In this initial step, a classic distance vector exchange routine 
takes place (as in most DV-hop type localization algorithms 
where all target nodes in the network get distances, in 
hops, to their reachable anchors i.e. at most K-hops away. This 
can be initiated by anchors nodes which broadcast beacons 
to be flooded throughout the network containing the anchor’s 
location with a hop-count value initialized to one. Each node 
maintains a table containing the coordinates (xi^pi) of anchor 
node i and the minimum hop distance to it hi, and exchanges 
such updates only with its immediate neighbours. 

B. Step 2: Multihop connection probability density function 

The conditional probability density function p{d\k) of inter¬ 
node distance d given the minimum hop count /c, is to the best 
of our knowledge the necessary ingredient which differentiates 
our algorithm from all others. From Bayes theorem we have 
that p{d\k) = p{k\d)p{d)/p{k). We may therefore calculate 
p{d\k) indirectly through p{k\d)p{d) and p{k). Since p{k) is 
independent of the inter-node distance, it does not affect the 
maximization calculation and so will be ignored later on. 

The conditional probability of the target-to-anchor hop 
count being k, given that they are a distance d apart is 
given by p{k\d) and can be approximately formulated in 
closed form as given by Gl), P^ . Similarly, the probability 
distribution of the distance between two random points p{d) 
can be formulated as in |T9| . However, the probability density 
function in (Tt) consists of multiple integrals and assumes 
very dense networks, making it difficult to calculate, even 
more so by the sensor’s weak processing unit. Therefore, 
in order to make progress we propose here a method for 
producing the probability density function p(k\d)p{d) and then 
flooding the network with the required fitting parameters. First, 
Monte Carlo simulations are employed to calculate the discrete 
values of p{k\d)p{d), which are then fitted into a continuous 
functions. 

1) Monte Carlo simulation: Let K < N and D ^ 1 
be the maximum allowed hop count and maximum allowed 
Euclidean distance between a target and anchor node respec¬ 
tively. Roughly, D should be at least larger than Kdmax 
or K/jS. Partition the disk of radius D into L concentric 
shells (like an archery target) of widths 26 = DjL indexed 
by I = 1,2...!/ such that each shell has central radius of 
di = 6{2l^l). We want to calculate the value of p{k\di)p{di) 
for each k = 1,2, ... if and each di with I = 1,2...!/. 
The Monte Carlo simulation can now be performed at one or 
more sensor nodes (preferably ones with significant processing 
power) or a central server e.g. the gateway as follows: 

1) Generate random coordinates with intensity p = N/A 
inside some predefined region. If the total number of 
nodes N and the WSN deployment region shape and 
area A are known then this is easy. If the WSN region 
shape is unknown, then a large square region can be 
used - it is shown that this does not affect the results 
significantly (see Fig. [^. If the density of nodes is 


unknown, then it can be estimated as described below 
in subsection IIII-B3I 

2) Generate communication links between nodes based 
on their mutual distances and appropriate connection 
probability function H (e.g. or ^ or other). 

3) For each k G [1,^] and I G calculate 

the discretized cumulative probability p{k\di)p{di) « 
P{k\{di - 6 < d < di + 6))P{di -6 <d<diP 6)/26. 

4) Repeat steps (1) - (3) above several times in a 
Monte Carlo fashion in order to refine the estimated 
p{k\di)p{di). 

2) Fitting: For each /c G [1, if], fit the discrete probability 
distribution p{k\d)p{d) to the following function of d 

pik\d)p{d) = eM-Mk)id - B{k)f + C{k)), (4) 

where A{k), B(k), C(k) are functions of k, e.g. polynomials 
of degree p ^ 1. In the simulations that follow we use 
p = 4. The Gaussianity of 0 was inspired by the extensive 
simulations results and analysis presented in pQ| . 

3) Density estimation: This section describes a simple 
method of estimating node density p, a prerequisite for per¬ 
forming the Monte Carlo simulations and building the said 
distributions. In a uniformly distributed network, the node 
density is defined as p = N/A, where A is the WSN 
deployment region area, thus it can be approximated by the 
average number of one-hop neighbours Ne divided by it’s 
average communication area A^. The former can be simply 
determined by each sensor counting the number of 1-hop 
neighbours from each node while the latter can be calculated 

p27T pOO 

Ae= / rH{r)drdO^ (5) 

Jo Jo 

where H{r) is the pair connectedness function of two nodes 
whose relative distance is r (see and The estimated 
density pe is then sent to the nearest processor which performs 
an average to obtain a refined estimate. 

C. Step 3: Maximum likelihood based Multihop localization 

This subsection describes our proposed kHopLoc algorithm 
utilizing maximum likelihood methods. First we introduce the 
likelihood function and then we address a method to maximize 
it. Consider a target node X and let hi be the hop count 
measure along the shortest path from X to anchor node i G 
[1, M]. The likelihood function of X at the yet undetermined 
coordinate (x, y) is defined as 

L{x,y)=p{x,y\hi,h 2 ,...hM)- (6) 

It follows that the best estimate of the true location of node 
X is the value that maximizes it’s likelihood 

{x*,y*) = axgmax L{x,y). 

{^,y) ^ ’ 

Assuming that the probability density functions p(x, y\hi) and 
p{x^y\hj) are mutually independent for i j, equation ([^ can 





(a) Isotropic network (b) Anisotropic network 

Fig. 2. Illustration of fixed anchor node locations. 


be written as 


M M 

L{x,y) = W_p{x,y\hi) = Y\‘^'^p{di\hi), 

i=l i=l 


( 8 ) 


where di = y/{x — Xi)^ {y — ViY Vi) the coor¬ 

dinates of anchor node i. Invoking Bayes theorem p{di\hi) = 
p{hi\di)p{di)/p{hi) yields 




i=l 


P{hi) 


(9) 


Substituting equation <0 back into equation Q , we obtain 


Y\tiPihi\di)-Y]::up{di) 

YltMhi) 


M 


(x*, ^*) = argmax 

{x,y) 


( 10 ) 


Because p{hi) is independent of {x^y), the product in the 
denominator can be eliminated from equation ( p^ such that 

M 

= argmaxTTp(/ii|di)p((ii). (11) 

{x,y) ZA 

Finally, substituting our fitted equation <0 into equation 
gives the result 


M 


(x*,^*) = argmaxTTexp(-A(/ii)(di - B{hi)f‘ + C{hi)) 

t=i 


M 


= argminy^ A{hi){di - B{hi)A. 
7^1 


( 12 ) 


The right hand side of equation (T?) can be easily calculated 
using gradient descent method or Newton method for example 
and can be performed by each target node independently. 


IV. Performance evaluation and analysis 

In this section, we demonstrate using computer simulations 
the effectiveness of kHopLoc, and also compare to the original 
DV-hop |[8| algorithm and a recent improved variant employing 
an approximate shortest path between nodes (H) which we 
will refer to as the ASP algorithm for short. The ASP algo¬ 
rithm is also a multihop range-free localization algorithm that 
tolerates network anisotropy with a small number of anchors. 
A detoured path detection is proposed which measures the 
deviation in the hop count between the direct and shortest 
paths of a node pair. A novel distance estimation method 
is introduced to approximate the shortest path based on the 



Fig. 3. Average localization error of DV-hop and kHopLoc in isotropic 
networks with random target node locations, a) N = 300 with random anchor 
node locations, h) N = 300 with fixed anchor node locations as in Fig.[^a). 
c) M = 13 fixed anchor node locations and N G [200, 700]. 


path deviation and to estimate their Euclidean distance by 
taking into account the extent of the detour of the approximate 
shortest path. 

To evaluate the performance of the kHopLoc in isotropic 
networks and anisotropic networks, we deploy sensor nodes 
randomly in a 10 x 10 square-shaped region (such that A = 
100) and in a 10 x 10 C-shaped of width 2 (such that A = 52), 
where communication probabilities in isotropic networks and 
anisotropic networks are assumed to follow Raleigh fading 
communication model and QUDG communication model, 
respectively. In Rayleigh fading communication model, we 
assume the path loss exponent to be = 2 and parameter 
/3 = 1 in formula In the QUDG communication model, 
we assume the maximum connection distance dmax = and 
the degree of irregularity DO I = 1.5. Moreover, to study the 
characteristic localization errors we simulate results in both 
isotropic and anisotropic WSN deployment regions with vari¬ 
able a) anchor node positions (fixed and random), b) number 
of anchor nodes, and c) node densities. Communication and 
computation overhead costs are also compared and discussed. 

A. Localization error 

Fig. depicts the average localization errors of DV-hop 
algorithm and kHopLoc by varying the number of anchor 
nodes and total sensor nodes in a square-shaped region under 
Rayleigh fading communication model. In all cases, the error 
based on kHopLoc is significantly smaller than that of the DV- 
hop algorithm with performance gains ranging between 20 and 
40%. Generally, the average localization error decreases with 
number of anchor nodes. Fixed anchor nodes (at strategic lo¬ 
cations as in Fig. |^a)) in general provides better performance 
than randomizing anchor nodes, especially when the number 
of anchor nodes is small. This is clearly due to the controlled 
avoidance of overlaps between anchor nodes which can lead 









































Fig. 4. Example isotropic WSN topologies and localization errors shown 
as disks. In all cases there are M = 13 fixed anchor nodes as in Fig [^a). 
Panels (a) and (c) use kHopLoc and N = 200, 700 sensor nodes. Panels (b) 
and (d) use DV-hop algorithm and N = 200, 700 sensor nodes. It is clear 
that kHopLoc has smaller localization errors, particularly near the boundary. 


to duplicate mutual information. 

Fig. 1^ illustrates two example network topologies (a sparse 
regime at = 200, and a dense regime at = 700), and 
highlights the localization error of each node. There are exactly 
M = 13 fixed anchor nodes in all cases. The green lines 
connecting the nodes and blue circles describe communication 
links and localization errors in which the radius are propor¬ 
tional to the localization errors. In kHopLoc (panels (a) and 
(c)), errors of nodes having few links tend to be big, thus the 
average error decreases in the dense regime of A^ = 700. On 
the other hand, in the DV-hop algorithm, localization errors of 
nodes near the border tend to be significantly larger. 

Fig. depicts the average localization error of the DV-hop 
algorithm |[^, the ASP algorithm II3)> and kHopLoc, for 
varying number of anchor nodes and total sensor nodes in 
C-shaped anisotropic region under the QUDG communication 
model. In all cases, the error based on kHopLoc is smaller 
than the other algorithms. Notice that there are two result 
curves for kHopLoc: when the density and shape of the 
deployment region are known (green curve), and when both 
the density and shape of the region are unknown (red curve). 
The performance of the case when the region shape is known 
is better than the other case since the Monte Carlo simulation 
phase of the algorithm over a known shape region produces 
more precise distributions and thus results for the multihop 
connection probability than when the region is assumed to 
be a square. Significantly, it is worth noting that the average 
localization error due to kHopLoc continues to decrease for 
the other two algorithms seems to saturate after 300 nodes. 
This demonstrates the benefits of using the full statistical hop- 
distribution generated during the first step of our algorithm 
(see section rather than just first order statistics such as 

the mean one-hop Euclidean distance. 

Fig. illustrates network topologies (a sparse regime at 
N = 200, and a dense regime at A^ = 700), and highlights the 



Fig. 5. Average localization error of DV-hop, ASP, and kHopLoc in 
anisotropic networks with random target node locations, a) N = 300 with 
random anchor node locations, b) M = 14 fixed anchor node locations (as 
in Fig.j^b)) and N e [200, 700]. 

localization errors of each node. There are exactly 14 fixed 
anchor nodes in the anisotropic C-shaped WSN deployment 
regions as shown in Fig. [^b). The green lines connecting 
the nodes and blue circles describe communication links and 
localization errors in which the radius are proportional to the 
localization errors. Similar to isotropic network, in kHopLoc, 
the main source of localization errors is due to nodes having 
few one-hop links. The reason for this is that nodes with fewer 
links tend to require a larger than average hop-count to reach 
anchors, thus making it difficult to estimate these node loca¬ 
tions accurately. As can be seen, this improves significantly 
in the dense regime. On the other hand, localization errors 
in the other two algorithms seem not to improve with node 
density. The reason for this is that these algorithms suffer from 
inaccurate inter-node distance estimations. 

B. Overhead Analysis 

This section discusses the communication and computa¬ 
tional costs of the three localizations algorithms under inves¬ 
tigation: DV-hop, ASP, and kHopLoc. 

1) Computational cost: For calculating location of one 
target node, DV-hop algorithm costs 0{M^) if using normal 
matrix multiplication. This cost comes from the matrix oper¬ 
ations necessary to estimated the least square error Q. The 
ASP algorithm costs consist of additional computation costs 
of 0{A) (where A is the number of tagged partitions for 
calculating Riemann sum) for calculating node density and 
M{M — l)/2 • 0(e“^/^) for estimating distances between 
the target node and M anchor nodes, where M{M — l)/2 is 
number of compound shortest paths, and 0(e“^/^) (in which 
e is the upper bound of the norm of the gradient) is the 
number of iterations of Newton method pT| fo r calculating the 
optimal central angles of the virtual holes ||l3|. Consequently, 
the total computation complexity of the ASP algorithm is 
0(M8 +A + .e-3/2). 

On the other hand, the cost of kHopLoc consists of Monte 
Carlo simulation cost and the MLE localization cost. The 
former computational cost requires generating the random 
network, finding the shortest k-hop paths and fitting. Thus, 
the cost comes up to 0{I • K • L • e“^/^), where / 

is number of Monte Carlo iterations and 0{N^) is the cost 
for calculating all-pairs shortest paths if using for example 

































(e) DV-Hop, 700 nodes 



(b) DV-Hop, 200 nodes 



(d) kHopLoc, 700 nodes 



(f) ASP, 700 nodes 


Fig. 6. Example anisotropic WSN topologies and localization errors shown 
as disks. In all cases there are M = 14 fixed anchor nodes as in Fig[^b). 


the Floyd Warshall algorithm and 0{K • L • is 

the fitting cost (see Q) if using Newton method. However, 
this computation can be done just once at a central node (or 
some server e.g. at the gateway) and then flooded through the 
network thus incurring an additional communication cost of 
0{N). Or otherwise, it is done before deploying the sensor 
nodes and then derived parameters (A{k), B{k)) are included 
into each sensor nodes, thus no additional communication cost 
occurs. 

The MLE cost is 0{M • where is the 

number of iteration of Newton method and 0{M) is cost 
for calculating the value of function in each iteration. 
Therefore the total cost amounts to 0(/ • + M • 

Obviously, the MLE localization cost is smaller than that of 
the ASP algorithm. 

2) Communication cost: The communication costs of the 
DV-hop algorithm and the ASP algorithm are bounded by 
20 {M{N — M)), where M and N — M are number of 
anchor and target nodes respectively. This is because these 
algorithms perform flooding twice - first for the minimum hop 
count estimation, and second to broadcast the average one-hop 
distance. 

The communication cost of kHopLoc however is mainly due 
to the initial hop count calculation giving 0{M{N — M)) (i.e. 
similar to |T^). When the density p is unknown, additional 
communication costs of 20(N) can be incurred, in which 
0{N) occurs when the nodes pass theirs 1-hop neighbour 


number to the central node. After the central processor runs 
the Monte Carlo simulations and fits the said distributions, 
it then broadcast the results (parameters A{k)^B{k)) to all 
nodes, thus costing another 0{N). The latter may lead to 
bandwidth issues in terms of the amount and resolution of 
the feedback information being fiooded, thus suggesting a 
trade-off between communication overheads and fitting and 
localization accuracy. 

V. Conclusions 

We have proposed a maximum likelihood based multihop 
localization algorithm called kHopLoc for use in wireless 
sensor networks (WSNs). The main advantage of the algorithm 
is the use of a Monte Carlo initial training phase to generate 
the multihop connection probability density functions. These 
are then used to build likelihood functions whose maxima 
estimate each target node location. Since the algorithm uses 
full statistical information for the multihop connection proba¬ 
bilities, localization results are significantly (about 20 — 40%) 
more accurate for both in isotropic and anisotropic networks. 
We have validated these results through computer simulations 
and discussed how and why some localization errors appear. 
Einally, we have discussed the communication and compu¬ 
tational costs of kHopLoc compared to conventional ones. 
Moreover, like most range-free algorithms, kHopLoc can be 
used in conjunction with GPS and/or range-based localization 
schemes to improve performance and energy consumption in 
WSNs. In the future, we aim to address the outstanding issue 
of localization in WSNs with non-uniform node deployments. 
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